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FOREWORD 


The  Global  Positioning  System  (GPS)  is  being  developed  primarily  for  the  purpose  of  providing  a  user 
with  accurate  three-dimensional  navigation  fixes  anywhere  near  the  surface  of  the  earth.  This  system  has  many 
unique  characteristics  that  make  it  useful  for  other  applications  as  well.  One  of  these  is  the  capability  of 
providing  the  user  with  orientation  information  if  he  is  able  to  receive  signals  from  several  antennas 
simultaneously. 

This  report  describes  a  preliminary  investigation  of  the  potential  for  simultaneous,  real  time,  navigation 
and  orientation  solutions  using  time-multiplexed  GPS  receiver  technology  such  as  the  Texas  Instruments  (TI) 
GEOSTAR.  Computer-generated  data  was  used  to  simulate  the  gross  characteristics  of  the  GEOSTAR.  This 
data  was  then  processed  with  a  least-squares  algorithm  to  provide  the  desired  solutions. 

Support  for  this  simulation  was  provided  by  the  Eastern  Space  and  Missile  Center,  Patrick  Air  Force 
Base,  Cocoa  Beach,  Florida. 

This  report  was  reviewed  and  approved  by  R.  W.  Hill,  Head,  Space  Flight  Sciences  Branch;  and  Carlton 
W.  Duke,  Jr.,  Head,  Space  and  Surface  Systems  Division. 
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IMPLEMENTATION 

Three-dimensional  orientation  requires  at  least  three  independent  antennas  to  define  a  plane.  Four  non- 
coplanar  antennas  could  be  used  to  define  two  planes  and  provide  redundancy.  Orientation  solutions  of  the 
planes  containing  the  antennas  are  related  to  the  vehicle  and,  therefore,  allow  the  orientation  of  the  vehicle  to 
be  determined.  The  vehicle  position  can  be  obtained  from  pseudoranges  in  a  manner  similar  to  the  conven¬ 
tional  navigation  so'ution  if  four  GPS  satellites  are  in  view1.  Four  satellites  will  also  add  strength  to  the 
orientation  solution.  For  these  reasons,  and  algorithm  simplicity,  this  simulation  requires  that  four  satellites 
always  be  in  view.  Operational  versions  of  this  software  will  need  to  be  able  to  handle  cases  where  fewer  than 
four  satellites  are  available,  but  this  complication  will  be  deferred  pending  the  results  of  the  current  simulation. 

The  TI  GEOSTAR2  receiver  has  eight  software  trackers  available.  These  are  intended  to  be  used  to  track 
four  satellites  on  both  L  band  frequencies  simultaneously.  The  two  frequency  measurements  provide  the  means 
to  correct  for  the  propagation  effects  due  to  ionospheric  refraction.  If  the  receiver  were  operated  on  a  single 
frequency,  there  would,  in  principle,  be  enough  trackers  available  to  allow  four  satellites  to  be  tracked  by  two 
antennas.  Assuming  that  the  antennas  are  not  separated  by  more  than  a  few  tens  of  meters,  the  loss  of  L2  would 
not  affect  the  orientation  solution  appreciably.  However,  the  position  solution  would  be  adversely  affected;  for 
this  reason  alone,  two-frequency  operations  will  probably  be  required. 

The  desirability  of  two-frequency  observations  from  four  satellites  implies  that  a  complete  GEOSTAR 
receiver  needs  to  be  dedicated  to  each  antenna.  Since  at  least  three  antennas  are  required,  three  GEOSTARs 
feeding  a  single  processing  computer  will  be  necessary.  This  is  probably  an  extravagant  use  of  hardware,  but  it 
may  be  necessary  to  realize  the  precision  desired. 

An  alternative  might  take  advantage  of  the  similarity  of  the  signal  dynamics  received  from  each  antenna. 
[  he  GEOSTAR  has  a  fundamental  clocking  interval  (T)  of  20  ms  duration.  All  receiver  operations  are  some 
integer  fraction  or  multiple  of  T.  Typically,  the  GEOSTAR  dwells  for  T/2  ms  on  each  satellite  and  T/4  on  each 
frequency  of  a  particular  satellite.  Thus,  it  has  completed  an  observation  cycle  (two  frequencies  and  four 
satellites)  after  2T  ms'.  If  the  antennas  were  connected  to  a  single  GEOSTAR  via  an  RF  switch  that  operated  in 
synchronism  with  the  GEOSTAR  cycle,  then  a  new  antenna  could  be  switched  on  line  every  2T  ms.  Since  the 
signal  dynamics  will  be  similar  from  each  antenna,  reacquisition  of  the  same  four  satellites  via  the  new  antenna 
using  the  GEOSTAR  receiver  software  should  not  be  difficult.  Thus,  it  might  be  possible  to  keep  three  or  four 
auxiliary  sets  of  tracking  loops  (one  per  antenna)  running  in  the  navigation  processor,  each  set  being  updated 
h>  the  receiver  processor  or  software.  Update  intervals  of  the  trackers  in  the  GEOSTAR  navigation  processor 
would  then  be  2NT  ms,  where  N  is  the  number  of  antennas  being  utilized. 

Several  variations  of  this  scheme  may  be  possible.  The  trackers  in  the  navigation  processor  might  be 
duplicates  of  the  receiver  processor  trackers,  except  that  they  would  be  updated  less  frequently.  They  could  also 
be  difference  trackers;  again  taking  advantage  of  the  similarity  of  signal  dynamics.  In  this  case,  one  antenna 
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would  be  the  reference  and  its  tracker  would  follow  the  full-signal  dynamics.  The  other  trackers  would  just 
track  the  difference  between  the  dynamics  received  by  its  antenna  and  the  reference.  These  difference  trackers 
would  not  need  frequent  updates  and  so  the  more  critical  reference  trackers  could  be  updated  in  alternate 
2T-ms  intervals.  The  antenna  switching  sequence  might  be  ABACADAB  .  . .  where  antenna  A  is  the  designated 
reference.  In  this  operating  mode,  the  GEOSTAR  navigation  processor  would  be  busy  tracking  signals  and. 
consequently,  an  outboard  computer  would  be  required  to  perform  the  real-time  position  and  orientation 
solutions. 

In  either  mode  of  operation,  each  tracker  will  be  updated  at  a  regular  but  different  interval.  Consequently, 
it  w  ill  be  necessary  to  propagate  the  current  state  of  each  tracker  to  a  common  future  time.  The  error  involv  ed 
in  this  operation  is  small  in  low-dynamic  applications,  and  it  has  the  important  benefit  of  providing  simultane¬ 
ous  observations  from  all  trackers  at  each  requested  observation  time.  The  simulation  uses  this  property  and. 
therefore,  requires  that  all  observations  occur  simultaneously  at  a  time  tag  obtained  from  the  local  clock. 

The  broadcast  ephemeris  and  satellite  clock  correction  from  each  satellite  are  used  to  compute  the  earth- 
fixed  coordinates.  These  ephemeris  parameters  are  assumed  by  the  simulation  to  be  correct.  There  are  no 
compensating  parameters  in  the  state  vector.  However,  the  local  clock  offset  is  unknown  and  is  a  parameter  in 
the  state  vector.  Also  unknown  are  antenna  biases.  These  states  model  the  differences  in  cable  length  between 
the  various  antennas  and  the  receiver.  These  differential  cable  delays  need  not  be  calibrated  as  long  as  they  can 
be  modeled  by  a  linear  polynomial  of  the  form 

p(t)  =  an  +  a,(t  -  t„) 

However,  it  is  assumed  that  there  are  no  unexpected  discontinuities  that  might  influence  the  coefficients  an  and 
ai  during  the  period  of  observation.  There  is  the  possibility  that  certain  other  effects  (e.g.,  interchannel  biases  in 
multihardware-channel  receivers)  may  behave  like  antenna  biases  over  the  short  term,  but  may  not  be  properly 
modeled  by  p(t)  over  the  long  term.  In  the  GEOSTAR,  there  is  only  one  hardware  channel  and  whenever  a  new 
satellite  is  acquired,  it  is  presumed  that  the  phac''  observations  from  the  new  satellite  are  consistent  w'ith  the  old 
and  not  randomized  by  an  arbitrary  phase  offset.  If  experience  demonstrates  that  this  condition  cannot  be  met, 
then  a  new  antenna  bias  will  have  to  be  defined  whenever  a  different  satellite  is  selected. 


CONDITIONS  OF  THE  SIMULATION 

I  he  observation  error  has  been  selected  on  the  basis  of  the  information  in  the  TI  Design  Analysis  Report 
III1  and  two-frequenev  operation.  The  GEOSTAR  observation  errors  follow: 

Pseudorange:  L[  =  0.687  m  RSS 
L;  =  0.768  m  RSS 

lonospherically  Corrected  Pseudorange  =  2.11  m 

Phase:  <P  Code  Tracking):  I.i  =  0.147  cm  RSS 

T;  =  0.199  cm  RSS 

lonospherically  Corrected  Phase  =  0.48  cm 
■"  6.48  oi  +  2.39  oi 
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All  simulated  observations  are  free  from  tropospheric  refraction  and  relativity  effects;  however,  flicker  noise 
has  been  added  to  model  the  expected  random  variations  due  to  the  local  frequency  standard.  The  simulation 
assumptions  follow; 

1.  Four  satellites  always  available 

2.  Observation  errors:  Pseudorange  a  —  2  m 

Phase  a  =  0.005  m 

3.  All  observations  simultaneous 

4.  Satellite  ephemeris  and  time  corrections  contain  no  error 

5.  Antenna  baseline  length  and  orientation  with  respect  to  vehicle  unknown 

6.  Antenna  biases  modeled  by  bias  and  drift 

7.  Tropospheric  refraction  and  relativity  neglected 

8.  Receiver  tracker  phase  does  not  randomize  when  switching  satellites 

The  GPS  satellite  constellation  used  in  the  simulation  was  the  18-satellite.  6-plane  system.  The  right 
ascension  of  the  ascending  nodes  of  the  constellation  are  separated  by  60  deg,  with  the  satellites  within  a  plane 
separated  by  120  deg.  The  phasing  of  satellites  in  adjacent  planes  is  offset  by  40  deg.'  The  orbital  period  is  the 
standard  1 1.966  hr,  eccentricity  w’as  selected  as  0.001,  and  the  inclination  of  each  plane  was  55  deg.  Figure  1 
shows  the  look  angles  from  the  receiver  to  the  five  satellites  in  view  during  the  simulation.  Of  these  five 
satellites,  7,  10,  15,  and  17  produced  the  best  Geometric  Dilution  of  Precision  (GDOP)  and  were  the  four  for 
which  synthetic  data  was  generated. 

ESTIMATOR 

The  data  processing  was  performed  by  a  least-squares  batch  estimator  using  both  the  pseudorange  and 
phase  data  from  all  antennas  at  each  observation  time.  Two  classes  of  states  are  employed  in  the  solution.  The 
first  class  includes  all  time-dependent  states:  the  vehicle  earth-fixed  coordinates  (x,  y,  and  z),  the  vehicle 
orientation  with  respect  to  the  local  vertical  (roll,  pitch,  and  yaw),  and  the  local  clock  offset.  These  states  are 
estimated  solely  from  the  data  obtained  at  a  single  observation  time  and  are  independent  of  previous  results  in 
all  respects  with  the  exception  of  the  correlations  that  exist  between  the  first  and  second  class  of  states. 

The  second  class  of  states  include  the  antenna  biases  that  are  assumed  to  be  independent  of  time.  This 
assumption  allows  all  previous  data  to  be  employed  at  each  observation  time  of  these  states.  The  accumulation 
of  data  is  accomplished  by  eliminating,  in  a  formal  elimination  procedure  at  each  observation  time,  the  first 
class  of  states.  The  results  of  the  elimination  are  then  saved  and  accumulated  through  the  entire  simulation.  The 
matrix  manipulations  leading  to  the  solution  are  described  in  Appendix  A. 

Since  the  vehicle  position  and  orientation  are  members  of  the  state  vector  and  the  data  is  collected  by 
several  antennas  with  positions  that  arc  not  states,  it  is  required  that  the  location  of  the  antennas  in  the  vehicle 
coordinate  system  be  accurately  known  beforehand.  The  assumption  requires  that  these  locations  can  be 
survcved  with  a  precision  better  than  the  phase  measurement  accuracy  so  that  they  do  not  contribute  to  the 
orientation  errors.  A  description  of  the  v;  rious  coordinate  transformations  that  are  required  is  presented  in 
Appendix  B. 
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AZIMUTH 
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5 

90 
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56.4 

7 
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80 

319.6 

30.4 

10 
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120 

55.0 

23.6 

15 

270 

40 

60.6 

71.7 

17 

330 

320 

187.1 

18.9 

FIGURE  1.  SATELLITE  SELECTION  FOR  SIMULATION 


SIMULATION 

Simulated  pseudorange  and  phase  data  were  generated  for  these  antennas  at  6-s  intervals  for  a  span  of  720 
s.  Vehicle  motion  was  at  a  constant  velocity  of  2.54  m/s  and  a  constant  heading  of  127.8  deg.  Sinusoidal  roll, 
pitch,  and  yaw  signals  of  different  periods  and  amplitudes  were  introduced  as  perturbations.  The  heading  was 
added  directly  to  the  yaw  signal.  The  local  clock  was  offset  from  GPS  time  by  50  ^s  with  a  drift  of  1  ns/s. 
Added  to  this  offset  was  the  flicker  noise  process  that  models  the  random  variations  of  the  local  crystal 
frequency  standard6.  Antenna  biases  were  also  introduced  on  all  antennas  except  the  reference.  The  biases  for 
the  reference  are  assumed  to  be  included  in  the  local  clock  offsets.  All  the  conditions  of  the  simulation  are 
summarized  below: 

Dynamic  Cases:  Time  Step  =  6  s;  Time  Span  =  720  s 

Baseline  System  Translation  and  Rotation:  Vhx  =  -0.61;  0,  =  -45 

Vhx  =  +  10.36  0:  =  -  5 
Vtv  =  +  20.73 
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True  States:  Longitude  =  282c  +  0.000018°  t 
Latitude  =  28°  -  0.000014°  t 
Height  =  0  m 

Roll  =  5  *  sin  ( 2tt  *  0.01)  t 

Pitch  =  10'  *  sin  (2rr  *  0.03)  t 

Yaw  =  8  *  sin  (2n  *  0.02)  t 

Clock  Bias  =  50  gts  +  I  ns  s  +  flicker 

Antenna  2  =  2.9919  m 

Antenna  3  =  5.9957  m 

Antenna  4  =  —2.9979  m 

Ricker:  t,  =  1:  r:  =  100;  r,  =  106;  r,  =  10  " 

The  results  from  two  cases  will  be  presented.  Both  use  all  the  same  conditions,  except  that  the  antenna 
baseline  lengths  are  different  by  a  factor  of  10.  In  both  cases,  the  antennas  are  located  in  a  coordinate  system 
offset  from  the  vehicle  by  the  vector  V.,  and  rotated  in  azimuth  by  0i  and  in  elevation  by  0:.  The  antenna 
locations  (meters)  in  this  rotated  system  are  as  follows: 


Case  1:  0,  0.  0  Case  2:  0,  0.  0 

2.  0.  2  20.  0.  20 

0.  2.  2  0,  20,  20 

RESULTS 

The  initial  conditions  given  the  estimator,  vs  the  true  conditions,  are  tabulated  in  Table  I .  The  given  con¬ 
ditions  were  offset  from  the  truth  by  2  deg  in  longitude  and  2  deg  in  latitude.  The  heading  or  yaw  angle  was  in 
error  by  37.8  deg.  These  offsets  were  included  to  demonstrate  the  rapid  convergence  to  the  true  solution  as  is 
illustrated  in  the  residual  plots.  Both  the  navigation  solution  and  the  orientation  solution  converge  within  the 
first  few  observation  intervals. 


TABLE  1.  INITIAL  CONDITIONS  OF  THE  SIMULATION 

True  Initial  C  onditions  Given  Initial  Conditions 


Longitude  (deg)  .  282  280 

Latitude  (deg)  .  28  26 

Height  (m)  .  o  0 

Roll  (deg)  .  o  () 

Pitch  (deg) .  o  0 

Yaw  (deg)  .  1 27.8  90 

Biases 

I  oca  I  Clock  lMs) .  50  0 

Antenna  2  (ns)  .  10  0 

Antenna  1  (ns)  .  20  0 

Antenna  4  (ns)  .  10  0 
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The  performance  of  the  estimator  is  evaluated  by  comparing  the  estimated  solution  for  each  state  with  the 
true  value  at  each  observation  time.  The  difference  (estimate-truth)  and  the  standard  deviations  of  various 
states  are  plotted.  The  average  values  of  the  residuals,  the  RMS  of  the  residuals  over  the  entire  720-s  span,  and 
the  standard  deviation  at  the  end  of  the  simulation  are  given  in  Table  2.  Plots  1-9  illustrate  the  residuals  from 
Case  1  and  Plots  1 1-19  are  residuals  from  Case  2.  The  plots  are  semilog  with  GPS  time  the  linear  horizontal 
axis  and  the  residuals  plotted  on  the  vertical  log  axis.  The  standard  deviations  at  each  time  are  plotted  as  heavy 
solid  lines  and.  in  most  cases,  bracket  the  residual  values  that  are  plotted  as  “o"  symbols. 

TABLE  2.  AVERAGE,  RMS,  AND  STANDARD  DEVIATION  OF  THE  RESIDUALS 
OVER  THE  ENTIRE  720-S  OBSERVATION  SPAN 


CASE  1:  3  ANTENNAS,  2-M  BASELINES 


State 

Average 

RMS 

S.D. 

j 

x  (m) 

0.096 

1.345 

1.315 

.  • 

v  (m) 

-  0.181 

2.052 

2.392 

z  (m) 

0.029 

1.324 

1.289 

roll  (rad) 

-  0.000262 

0.002004 

0.002063 

pitch  (rad) 

0.000336 

0.002474 

0.001934 

yaw  (rad) 

-  0.000114 

0.002374 

0.002671 

clock  (m) 

-  0.040 

1.354 

1 .533 

CASE  2:  3  ANTENNAS,  20-M  BASELINES 

. 

x  (m) 

0.090 

1.340 

1.314 

*  •  * 

y  (m) 

-  0.185 

2.052 

2.392 

z  (m) 

0.026 

1.321 

1.289 

roll  (rad) 

+0.000035 

0.000386 

0.000206 

pitch  (rad) 

0.000245 

0.000739 

0.000193 

yaw  (rad) 

+0.000076 

0.000610 

0.000267 

clock  (m) 

-  0.044 

1.354 

1.533 

If  the  estimator  has  done  its  job  properly,  the  residual  plots  of  all  Class  1  states  should  contain  no 
systematic  signals  or  trends.  If  the  estimator  has  modeled  the  physical  system  correctly  ,  the  residuals  should  be 
representative  of  the  noise  processes  that  are  inherent  in  the  measurement  and  reflected  in  the  states.  Class  2 
states,  the  antenna  biases,  are  correlated  Irorn  observation  to  observation  due  to  the  accumulation  process. 
Consequently,  the  Class  2  residuals  show  considerable  correlation  as  they  hunt  lor  the  zero  residual  level. 

Plots  |-3  and  I  1-13  show  the  nav  i  gat  nm  residual -  m  the  twn  eases  !  he  change  in  baseline  length  has  little 
et  led  upon  the  navigation  solution  or  upon  the  h  'cal  i  .<  k  k  ot  (set  (  Plots  and  I  ).  1  it  is  result  is  expected,  since 
nav  i  gat  ion  can  be  perl  or  rued  with  a  single  antenna  I  1  oi  .entat nm  ic'iduals  are  shown  m  Plots  4-6  and  14-16. 
I  he  improv  ement  w  it  h  the  longer  baseline  I  (  a  si  7 1  ■  ad  .  apparent  I  he  standaid  deviations  ol  the  three 

orientations  approach  0.002  rail  in  (  a  sc  I  and  niioo.''  ;.u!  n  (  asc  2  Mi  cmploimg  the  assumptions  ol  this 

I: 


NSWC  TR  83-377 


The  construction  of  this  A  matrix  is  pursued  in  the  next  section  of  this  appendix 


AAx  k  =  (□  -  C)k 

A'WAAx  k  —  A'W  (a  -  C)k  (A-2) 

The  matrix  W  is  a  weight  matrix  that  weighs  the  two  data  classes  appropriately.  This  expression  can  he 
simplified  by  letting 


and 


B  -  A'W  A 
E  =  D  -  C 


Next,  the  B  matrix  is  partitioned  to  allow  for  two  classes  of  states:  Class  I  is  represented  by  the  states  Ax,,  and 
Class  2  by  the  states  Axh.  The  partitioning  is  as  follows: 


BppAxr,  +  BPh  Axi,  —  EP 

(A-4) 

Bph  Ax  p  +  Bsp  Ax h  =  Eh 

(A-5) 

Class  I  states  are  eliminated  by  solving  equation  (A-4)  for  Ax  ,: 

Ax  ..  =  B,,'  (EP  -  B,,h  Axh) 


Substitution  ot  this  resiilt  into  equation  (A-5)  produces  the  matrix  equation 

<B  B1  K  Axi.  =  Ek-  Bp„  B..J,  Er  (A-6) 

I  his  equation  •>  independent  ol  the  Class  I  states  Ax...  It  equation  (A-6)  is  summed  over  all  observations  up  to 
the  ptesetit  it  carries  .ill  the  inhumation  acquired  about  the  Class  2  states. 


IT  ,  1  j  B  B  Bp,  Bo-  | 

»  ;  ,  -  I  H'  B  p.  I  1 

I  hese  accumulated  matrices  can  then  he  used  along  with  the  data  from  the  most  recent  observation  at  t.  to 
estimate  the  current  state  I  his  linal  form  is  as  tollows 
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APPENDIX  A 


BATCH  LEaSE  SQUARES  FORMULATION 


MATRIX  MANIPULATIONS 


The  observations  used  in  this  simulation  consist  of  (1)  pseudoranges  from  each  GPS  satellite  to  each 
antenna  and  (2)  phase  differences  between  each  antenna  and  the  reference  antenna  from  each  of  the  satellites.  If 
there  are  n.  satellites  and  n,  antennas,  the  number  of  observations  at  each  time  step  tk  is 

n.  (2n.,  -  1 ) 

At  each  time  step,  there  will  exist  an  estimate  of  each  observatio.  based  upon  previous  information.  If  the 
actual  observations  are  represented  by  the  column  vector  Qk  and  tne  estimate  of  the  observations  by  Ck.  then 
the  difference  provides  information  about  the  error  in  the  estimate.  The  parameters  of  interest  are  the  states  Xk. 
and  since  Ck  is  computed  from  estimates  of  these  states  at  each  time  step,  the  error  in  the  estimate  provides 
information  to  correct  the  states.  The  formulation  that  adjusts  the  states  to  match  the  observations  begins  by 
expanding  Ck  by  use  of  Taylor’s  Series  around  the  current  state  estimate  Xk.  The  general  expression  for  the 
expansion  follows: 


dC.  -  I  -  c.  -  .  I  anCk 

-  c.  -  —  -  -  (X.  X.)  *  -  (X,  -  xK )  +  .  .  .  +  —  —  (xk  -  xj 

ax.  2'  dx.'  n!  0Xk 

In  order  to  employ  linear  estimation  theory,  this  expression  must  be  approximated  by  truncating  the  series 
alter  the  term  involving  the  first  partial  derivative  I  hen.  C.  is  moved  to  the  left  hand  side  of  the  equality  and 
what  temains  is  the  linear  approximation 


(  . 

□  .  (.  t-  (x.  x.t  ( A-I ) 

I  >  n.a'Mx  o|  pa  Mia!  del  natives  C.  X.  1  is  generally  known  as  the  “A”  matrix: 

•(  ac  ac  ac 

dx  d\:  a.x,  dx,  k 

I  ■  ' ,i hser ipt  !  represents  the  number  of  states  in  the  solution  and  k  is  the  time  step.  The  difference  between  the 

one  and  estimated  state  will  be  designated  by  the  vector  Ax  k.  where 

Ax„  =  itk  -  xk 
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GPS  TIME  (s) 

PLOT  18.  ANTENNA  2  BIAS  RESIDUALS  (ESTIMATE-TRUTH); 
20-m  ANTENNA  BASELINES 


PLOT  19.  ANTENNA  3  BIAS  RESIDUALS  (ESTIMATE-TRUTH); 
20-m  ANTENNA  BASELINES 
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o  °o 

% 

'  <?  ° 

vJ 

o  o  _ 

U  U  _ 

_ _ 2 _ 

GPS  TIME  (s) 

PLOT  7.  LOCAL  CLOCK  OFFSET  RESIDUALS,  (ESTIMATE-TRUTH); 
2-in  ANTENNA  BASELINES 


GPS  TIME  (s) 

PLOT  8.  ANTENNA  2  BIAS  RESIDUAL,  (ESTIMATE-TRUTH); 
2-m  ANTENNA  BASELINES 


GPS  TIME  (s) 
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model  fits  the  cable  behavior.  Unpredictable  phase  variations  due  to  kinks  or  thermal  effects  of  sunlight  and 
shadow  are  not  modeled  and  must  be  examined  separately. 

Clearly,  a  larger  baseline  pn  '.ices  higher  precision  orientation  capability  for  a  given  measurement  accu¬ 
racy.  However,  as  noted  above,  longer  baselines  are  prone  to  introduce  unmodeled  effects  that  may  make  the 
potential  improvement  unrealizable.  This  simulation  demonstrates  that  the  three-antenna  array  nominally 
separated  by  20-m  approaches  0.2  mrad  resolution  in  the  three  angles  for  the  given  0.5-cm  phase  measurement 
accuracy.  As  expected,  the  2-m  array  is  a  factor  of  10  worse.  In  general,  the  angular  precision  can  be 
summarized  by  the  following  relationship: 


A  t/r 


where  A  i/t  =  angular  precision  (rad) 

A  1  =  phase  measurement  error  (m) 
b  =  baseline  length  (m) 

i It  =  angle  between  the  baseline  vector  and  the  vector  to  the  satellite 

Use  of  several  satellites  and  two  or  more  baselines  will  tend  to  minimize  the  effect  of  eset/r  in  this  equation. 
Thus,  A  l,  b  more  nearly  represents  the  standard  deviation  observed  for  roll,  pitch,  and  yaw  in  this  simulation. 


CONCLUSION 

All  troublesome  propagation  effects  have  been  neglected  in  this  simulation  and,  consequently,  the  results 
represent  a  best  case.  Under  these  optimistic  conditions,  the  simulation  shows  that  milliradian  angular  solutions 
are  feasible  using  the  GPS  and  a  suitable  receiver  such  as  the  T1  GEOSTAR.  Submilliradian  accuracies  are  also 
possible,  but  will  require  a  distance  between  antennas  of  about  20-m.  It  is  doubtful  that  additional  improve¬ 
ments  in  phase  measurement  accuracy  will  buy  a  proportional  amount  of  improvement  in  angular  resolution. 
Effort  should  instead  be  placed  toward  improving  the  models  of  propagation  effects,  working  to  remove  the 
mechanical  instabilities  of  the  antenna  baseline  on  a  moving  vehicle,  and  controlling  the  thermal  variations  in 
long  signal  cables.  Clever  mechanical  design  and  more  sophisticated  software  will  be  needed  to  make  0. 1-mrad 
orientations  routine. 
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simulation,  longer  baselines  would  continue  to  improve  the  angular  resolution.  However,  practical  considera¬ 
tions,  discussed  in  the  next  section,  would  work  in  the  opposite  direction. 

The  antenna  biases,  initially  unknown,  gradually  converge  toward  a  small  residual  as  data  is  accumulated. 
The  bias  residuals  appear  in  Plots  8-9  and  18-19.  As  discussed  previously,  the  knowledge  of  these  states 
improves  due  to  the  accumulation  of  information  over  the  span  of  observations.  The  effect  of  this  improvement 
is  reflected  in  the  orientation  residuals  and  can  be  seen  particularly  in  the  case  of  pitch  (Plot  15).  During  the  first 
2  min  of  observations,  the  pitch  residual  is  larger  than  it  is  later  when  the  antenna  bias  residuals  are  reduced.  It 
can  also  be  seen  in  Plot  5.  but  the  larger  residuals  make  the  effect  less  obvious.  In  both  cases,  the  pitch 
component  shows  the  effect  more  than  roll  or  yaw  because  the  correlation  coefficients  of  pitch  with  respect  to 
the  antenna  biases  have  the  same  sign.  The  correlation  coefficients  of  roll  and  yaw  with  respect  to  the  antenna 
biases  have  opposite  sign  and.  consequently,  tend  to  cancel  each  other  when  the  antenna  bias  residuals  are  in 
the  same  sense. 


DISCUSSION 

This  simulation  demonstrates  that  the  intrinsic  measurement  accuracy  of  the  GPS  signals  readily  meets  the 
20-m  three-dimensional  navigation  error  requirement  of  the  application.  However,  there  are  many  sources  of 
measurement  error,  not  included  in  the  simulation,  that  are  capable  of  biasing  the  result  and  its  effects  are 
difficult  to  detect. 

Accurate  prediction  of  the  satellite  trajectory  over  a  long  enough  period  of  time  with  sufficient  accuracy  to 
make  the  broadcast  ephemeris  suitable  for  high-precision  navigation  and  point  positioning  will  continue  to  be  a 
challenge.  Improvement  of  our  understanding  of  the  forces  acting  on  the  satellites  and  an  increase  in  the 
number  of  tracking  stations  may  combine  to  decrease  the  size  of  the  broadcast  ephemeris  error  in  the  future. 
Closely  related  to  the  ephemeris  error  is  the  prediction  of  satellite  time.  Satellite  clock  corrections  are  included 
in  the  broadcast  ephemeris  and  need  to  be  considered  as  part  of  the  ephemeris  prediction  problem. 

Ionospheric  refraction,  tropospheric  refraction,  and  relativity  are  biases  that  were  not  considered  in  this 
simulation.  Errors  here  will  particularly  effect  the  navigation  solution, since  they,  like  ephemeris  errors,  can 
introduce  time-varying  radial  signals  that  cannot  be  easily  separated  from  the  true  signal.  The  orientation 
solution  is  not  significantly  degraded  by  these  types  of  errors  as  long  as  the  phase  front  remains  intact  across  a 
dimension  comparable  to  the  longest  antenna  baseline  length.  However,  warping  of  the  phase  front  due  to 
propagation  effects  so  that  arrival  times  at  the  antennas  are  perturbed  would  be  of  serious  concern  . 


Anv  phenomena  that  differentially  influences  the  arrival  times  of  the  signals  from  the  several  antennas  at 
the  receiver  in  an  unexpected  manner  would  produce  errors  in  the  orientation  solution.  Any  movement  of  the 
antenna  phase  center  as  a  function  of  the  direction  to  a  satellite  must  be  measured  and  modeled.  The  distance 
between  antenna  phase  centers  (i.e..  the  respective  baselines)  and  the  orientation  of  the  antenna  array  with 
respect  to  the  vehicle  must  all  be  measured  and  remain  stable  to  an  accuracy  better  that  that  desired  for  the 
system  performance.  This  requirement  tends  to  favor  short  baselines  and  systems  where  all  antennas  are  rigidly 
mounted  on  a  single  structure  that  can  then  be  precisely  mounted  on  the  vehicle. 

Propagation  delays  of  signals  from  each  antenna  to  the  receiver  can  be  accounted  for  in  the  antenna 
biases,  if  the  effects  are  time-independent.  Calibration  of  these  signal  cables  will  not  be  necessary  as  tongas  the 
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R[»i*k  A  x i*k  B[iUk  A  x  hk  E,,k 
B,Ii.k  Ax,,k  +  (Bi,hk  +  Dk-i)  Axi,k  =  (Ehk  +  E\.|) 


( A-7) 


Solution  at  tk  proceeds  by  premultiplying  both  sides  of  the  unpartitioned  matrix  equation  by  B  The  result  is 
the  vector  Axk  as  desired. 


i 

C  ONSTRI  CTION  OF  THE  A  MATRIX 


As  discussed  above,  two  data  classes  are  used  to  produce  the  combined  position  and  orientation  solution. 
The  partial  derivatives  for  each  of  these  will  be  presented  independently.  Reference  to  the  diagram  in  Appendix 
B  (Figure  B-l)  will  be  helpful. 

Pseudorange 


The  pseudorange  observations  can  be  modeled  by  the  following  equation: 

Pkn  -  |T;  -  e,  |k  +  C  [ro  +  r,  +  r,  (tk  _  to)]  i  >  I  (A-8) 

where  subscript  j  =  the  satellite 

subscript  i  =  the  antenna 

The  two  vectors  e,  and  r,  (defined  in  Figure  B-l)  are  the  ECEF  coordinates  of  the  i  th  antenna  and  the 
ECF.F  coordinates  of  the  j  th  satellite,  respectively.  The  biases,  represented  by  the  taus  in  equation  (A-8). 
are  composites.  The  local  clock  bias  r„  and  the  reference  antenna  biases  n  and  f  i  are  grouped  together  to 
make  a  single  Class  I  state  The  Class  2  states  are  the  biases  and  drifts  associated  with  the  remaining 
antennas  plus  (u.  These  effects  arise  due  to  cable  variations,  etc.,  and  are  represented  by  r,  and  r,. 
respectively  (i  >  I).  The  sum  of  all  these  effects  is  the  total  delay  impressed  upon  the  signal  from  each 
antenna.  In  summary,  the  biases  are  as  follows: 

reference  antenna:  -  n,  +  r,  +  r,  (tk  -  t») 

all  other  antennas:  +  r,  +  r,  (tk  -  to) 

Equation  (A-8)  is  the  observation  equation  for  the  pseudorange  data  class.  It  represents  Cu  which  was 
defined  in  the  first  section  of  this  appendix.  In  order  to  form  the  A  matrix,  the  partial  derivative  of  pk|,  is 
required  with  respect  to  each  state  variable: 

e„i.  e„;,  c\,<.  7i,  7:.  71.  £»,  r . r„.,.  r . r„., 


The  partial  derivatives  with  respect  to  these  states  are  written  below,  with  the  subscript  m  indicating  the 
coordinate,  m  l.  2.  3. 
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Spki. 


,  ,  3cxi  3ev,  /  3e„ 

(rXJ  —  e„ )  -jr- —  +  (ry,  -  e>,)  +  (r„  -  e„) 


dy„ 


dym 


dyn 


dym 

9pk.ii 


I  r,  —  e,  |k 


3pk 


afo 

3pku 


3r, 


=  c  i  >  1 


dr', 


=  c(tk  -  to)  i  >  1 


The  result  from  Appendix  B,  equation  (B-l)  is  required  to  evaluate  the  partials  of  the  antenna  positions 
with  respect  to  the  states.  The  partial  of  equation  (B-l)  with  respect  to  coordinates  follows: 


9e, 


3eu 


Teu 


3T„ 


3e„ 


+ 


3Te: 


3e„r 


r„,  Tvh  b, 


Since  T,„  is  independent  of  e;.  (Figure  B-3),  the  above  matrix  equation  reduces  to 

3e,  3T,„ 


3e.„ 


3e„rr 


TUI  T,h  b, 


m  =  1,  2.  3 


The  evaluation  of  3T.,„  3e„„,  is  presented  in  Figure  A-l;  the  form  of  TC„.T,„.  and  T,h  are  given  in  Appendix  B. 


The  partial  of  equation  (B-l)  with  respect  to  the  orientation  angles  roll,  pitch,  and  yaw  follows: 


3e 


-  -  T,„ 

m 


3T„, 


3T,„ 


T,„  T,h  b, 


9*yrn  L  3  "y  m 

In  this  case.  3Teu  dym  -  0,  and  this  allows  the  above  matrix  equation  to  be  simplified  to 


3e, 


37m 


I 


=  T, 


3T,n 

37m 


T.b  b, 


m  =  1.  2.  3 


The  evaluation  of  3T„,  37  , m  is  presented  in  Figure  A-2.  All  these  pieces  can  be  put  together  to  form  the 
required  partial  derivatives. 

I  Phase  Difference 


The  observation  equation  for  the  phase  difference  data  class  can  be  obtained  by  subtracting  the  phase  predicted 
at  the  reference  antenna  from  the  phase  at  all  the  other  antennas: 

|  Apk„  -  p k pi  ~  Pkii  i  >  I  (A-9) 


The  expansion  of  the  above,  using  equation  (A-B).  is 
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APPENDIX  B 

COORDINATE  TRANSFORMATIONS 

The  GPS  broadcast  ephemeris  provides  the  user  with  the  satellite  position  in  the  earth-fixed  coordinate 
system'.  To  be  able  to  perform  the  least  squares  differential  correction  algorithm  described  in  this  report,  one 
must  be  able  to  calculate  the  receiver  antenna  positions  in  the  same  system.  The  locations  of  the  antennas  are 
assumed  to  be  surveyed  in  some  convenient  local  reference  frame.  The  relationship  of  this  frame  to  the 
earth-fixed  frame  is  the  subject  of  this  appendix. 

The  diagram  in  Figure  B-l  shows  two  antennas:  one  at  Bt  and  one  at  B,.  One  satellite  is  at  P  in  the  earth 
centered  earth-fixed  (ECEF)  coordinate  system.  The  vector  e,  represents  the  coordinates  of  the  i  th  antenna 
location,  the  vector  7,  represents  thej  th  satellite  location,  and  the  baseline  vectors,  connects  the  reference 
antenna  (I  =  1)  to  the  i  th  antenna.  These  relationships  are  described  mathematically  in  the  figure. 

The  observations  consist  of  pseudoranges  and  phase  differences.  The  pseudoranges  include  the  geometric 
ranges  h„  plus  the  antenna  biases  r,  and  the  local  clock  offset  r„.  The  phase  differences  include  the  segment  /„ 
plus  the  difference  in  antenna  offsets  r ,  -  r\.  These  observations,  obtained  from  three  or  more  antennas  and 
four  satellites,  are  required  to  produce  a  solution. 

The  antenna  locations  are  defined  in  a  reference  frame  independent  of  the  vehicle.  This  reference  frame  is 
the  bi,  b;.  b,  frame  drawn  in  Figure  B-2.  Each  antenna  has  coordinates  B,  (b,,.  b;,.  bi.)  in  this  baseline  system. 
The  vehicle  reference  frame  is  represented  by  the  v,.  v:,  Vi.  coordinates.  The  baseline  system  is  rigidly  attached 
to  the  vehicle  system,  but  is  translated  by  vh  and  rotated  by  the  angles  0i  and  0?.  The  vector  vb  and  the  angles  0, 
and  0:  are  assumed  to  be  known.  The  transformation  matrix  Tlh,  which  relates  any  point  in  the  baseline  system 
to  its  equivalent  in  the  vehicle  system,  is  presented  in  Figure  B-2. 

In  a  quiescent  state,  the  vehicle  system  and  the  local  vertical  system  would  be  coincident.  Flowever.  under 
dvnamic  conditions,  roll,  pitch,  and  vaw  are  permitted;  therefore,  the  vehicle  system  can  rotate  about  each  axis 
as  illustrated  in  Figure  B-3.  These  rotations  are  the  orientation  angles  of  the  vehicle  system  with  respect  to  the 
local  vertical  and  are  states  of  the  solution  vector.  The  transformation  matrix  Tm,  which  transforms  any  point 
in  the  vehicle  system  to  its  equivalent  in  the  local  vertical  system,  is  also  presented  in  Figure  B-3.  Note  that  the 
first  three  rows  of  the  fourth  columns  are  zero.  This  indicates  that  there  is  no  translation  and  so  the  origins 
always  coincide. 

A  third  transformation  relates  the  local  vertical  to  the  ECEF  coordinate  system.  Figure  B-4  shows  that  the 
local  vertical  system  is  right  handed  with  the  u(  axis  positive  in  the  radial  direction  and  u,  positive  eastward.  It 

A  .1  Van  Dierendonck.  S.S.  Russell.  E.R.  Kopit/kc.  and  M.  Birnbaum.  “The  GPS  Navigation  Message".  NAVIGA¬ 
TION.  Vol  25.  «2.  I97X 
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is  translated  from  the  origin  of  the  ECEF  system  by  the  vector  e„  and  rotated  by  X  and  /3.  The  e„  vector  joins  the 
origin  of  the  two  coordinate  systems  and  is  different  from  e,  in  Figure  B-l,  which  joins  the  origin  of  the  ECEF 
system  to  each  antenna.  The  angles  X  and  are,  respectively,  the  longitude  and  latitude  of  the  vehicle  origin 
and  can  be  computed  if  e„  is  know  n.  The  elements  of  the  vector  e„  are  state  elements,  and  along  with  roll,  pitch, 
yaw,  and  the  local  clock  bias,  complete  the  first  class  of  states  described  in  Appendix  A.  The  transformation 
matrix  Ttl„  which  relates  any  point  in  the  local  vertical  system  to  its  equivalent  in  the  FCEF  system,  is  shown  in 
Figure  B-4. 

These  three  transformations  together  transform  the  local  antenna  coordinates  b,  into  the  ECEF  coordi¬ 
nate  system  e,: 

e,  =  Tcu  T„,  T,h  b,  (B-l) 


z 


h  „  =  (r„  -  e,,)  x  +  (rM  -  ej  y  +  (r„  -  e„)  / 


s  ■=  (exl  -  evi)  x  +  (ex,  -  ex i)  v  +  (e„  -  e,,)  /  .  i  ^  1 

h;,  =  1,  +  h;  [ 
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FIGURE  B-4.  TRANSFORMATION  T«„:  LOCAL  VERTICAL  IN  ECEF  SYSTEM 
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